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Abstract. Marcus et al. (Marcus P, Ma H and Qiu S L 2002 J. Phys.: Condens. 
Matter 14 L525) claim that thermodynamic properties of materials under pressure 
must be computed using the Gibbs free energy G, rather than the internal energy 
E. Marcus et al. state that "The minima of G, but not of E, give the equilibrium 
structure; the second derivatives of G, but not of E, with respect to strains at the 
equilibrium structure give the equilibrium elastic constants." Both statements are 
incorrect. 
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Marcus et al. |T| presented an analysis of structural and elastic properties of solids 
subject to compression under athermal (T = 0) conditions. They claim that the Gibbs 
free energy G must be used at finite pressure P to find the equilibrium structure rather 
than the internal energy E. In particular, they consider the epitaxial Bain path (EBP) 
that relates body-centered cubic (bcc), body-centered tetragonal (bet), and face-centered 
cubic (fee) structures for Fe at 100 GPa (see also ref. [2J). They show that while 
G(c/a;p) along the EBP (with a and c the two independent lattice parameters) yields 
a minimum at c/a = 1 (the bcc structure), the minimum of E(c/a;p) along the EBP 
is displaced, at c/a = 0.95. This is a result of misusing elementary thermodynamics 

There is a minimum principle for the internal energy E at constant entropy S and 
volume V, for the enthalpy H = E + PV at constant S and P, for the Helmholtz free 
energy A = E — TS at constant temperature T and V, and for the Gibbs free energy 
G = E — TS + PV at constant P and T. There is no minimum principle for E (or 
H, which are equivalent at T = 0, considered by Marcus et al.) at constant P. The 
correct analysis is shown in Fig. 3 of Stixrude et al. [I], where E versus c/a is shown at 
constant V for the EBP for Fe, and the extrema are at the bcc structure, as expected. 
There is nothing wrong with minimizing G at constant P and T, but the exact same 
results will be obtained for minimizing E at constant V and S. The resulting pressure 
P can always be obtained from E, since P = — |^ \$- At T = the constraint dS = 
is trivial, since S = at T = 0. 

Marcus et al. further assert that proper elastic constants must be calculated from 
second Eulerian strain (e) derivatives of G (defined as Gy in pQ) rather than the 
internal energy E (c%). Calculating the shear elastic constants d = (c n — ci 2 )/2 
and C44 for fecc-iron using both G and E (c' and C44) they find a shear instability 
at 150 GPa using G (c') but not for E (c'), implying that previous computational 
estimates of elastic constants at pressure are incorrect, and that pressure corrections 
need to be applied. To the contrary, computation of elastic constants as prescribed 
by Marcus et al. gives incorrect results. Marcus et al. ignored the fact that the 
pressure and shear stresses vary as a function of strain. They did not obtain any 
thermodynamically valid second derivatives by their finite difference procedure, in 
which they computed ^ (Ei(eij) + PoV^Sy) — E — PqVo); they did not even obtain the 

derivatives y^^p- \e k 11 which are not in any case elastic constants, since the pressure 
and shear stresses vary with their deformation e^. This problem remains even if the 
second order coefficients of a fit for a polynomial expansion in y- (E + PoV) is obtained. 
The high pressure elastic constants computed using their procedures (21 El El 13 E] are 
incorrect, although the resulting errors may be small in some cases. 

It is important to use the appropriate thermodynamic function for the appropriate 
conditions, and as any student of elementary thermodynamics knows, it is also important 
to keep track of what is being held constant for a given partial derivative. Elastic 
constants can be defined by various ways: (1) from the equations of motion (i.e. 
sound velocities), (2) as derivatives of stress with respect to strain, or (3) as second 
derivatives of the internal (giving the adiabatic elastic constants) or Helmholtz free 
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energy (giving the isothermal elastic constants) with respect to strain, holding the other 
strains constant Pj. Different constants can also be derived depending on the use of 
finite or infinitesimal strain parameters. All of these definitions are equivalent at zero 
pressure, but differ under applied stress. The different definitions of elastic constants 
under applied stress remain sources of confusion ^Oj- Under no conditions, however, 
is an elastic constant tensor properly defined as derivatives of H or G with respect to 
strains, holding the other strains constant. 

Kamb ^1] even comes to the conclusion that it is not possible usefully to associate 
a Gibbs free energy with a non-hydro statically stressed solid; this is similarly stated by 
Wallace ^2]- The definition of elasticity on the basis of G is not well-founded. A 
problem arises, for example, when considering phase equilibria of a fluid in contact with 
a crystal surface; the chemical potential of components in the fluid in equilibrium with 
the solid vary according to the crystal face for a stressed solid, indicating there is no 
unique definition of the Gibbs free energy for a stresses solid [11 . 

Marcus et al. state that the elastic constants in Refs. f3J El D2] are incorrect and 
require pressure corrections. The elastic constants presented in Refs. [T3J E3 EEB \U\ 
are the elastic constants for wave propagation, and these are most easily measured, and 
are important in seismology and other applications. 

For isotropic initial stress the elastic constants for acoustic wave propagation and 
stress-strain coefficients are equivalent (see section 5 in jH]). We will now illustrate that 
the expression of strain-energy density can give the same elastic constants as the stress- 
strain relations for volume conserving strains for a reference state with isotropic applied 
stress. We use the fourth rank tensor notation from [5] for the elastic constants c^i 
(the stress-strain coefficients). Consider the expression for strain-energy density from 
Barron and Klein 

V E a + ~ (cijki - t,P (2$ij$ki - $u$jk - $ji$ik)) £ij£kh (!) 



V 2 \ J 2 

where 5ik is the Kronecker delta. Evaluating this expression, for example, for C1313 
(corresponding to C55 = C44 in Voigt notation, as used in ref. PQ), with the strain 

/ d \ 

e (d) = 0, (2) 
v d i 

with d the strain amplitude, one does indeed obtain a pressure correction term 

-y-2 [c 1313 + -pj d\ (3) 

This is the strain used in ref jHj, and their elastic constants should consequently be 
corrected for pressure to obtain wave propagation velocity. 

However, by choosing specific volume conserving strains, such as those given in 
[Tol IT7j . corrections can be avoided. For the C1313 example we apply the monoclinic 
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eq. [U becomes 
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For hexagonal and tetragonal systems C1322 
right hand side is zero. Expanding this into a series of d yields 

5 2 (2c 131 3) + 0[n 



2c 



d 3 



2 — (1-^2)2 ' — 1322- 1 _ d2 . 

and the final term in the sum on the 
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without any pressure correction term. This is also true for the other strains given in 

Halm. 

Single crystal elasticity is difficult to measure in high pressure experiments, 
especially for opaque materials such as metals which are discussed in [I]. However, 
advances in optical spectroscopy have made it possible to measure the Raman active 
phonon mode in hep metals. This optical mode and the shear elastic constant C1313 
can be viewed as properties of the same, nearly continuous phonon branch in an 
extended Brillouin zone scheme, and can be related by a simple force constant model 
of phonon dispersion [THj. Results obtained for Fe and Re compare very favorably 
with experimental estimates ^01 EDj over a wide pressure range, corroborating that no 
pressure correction need be applied. Computed finite temperature elastic constants for 
Ta [2T] also agree well with sound velocities obtained under shock conditions along the 
Hugoniot. 

To conclude, the internal energy E in conjunction with its volume V and strain 
derivatives defines the thermodynamics and elasticity of a material completely, even 
under applied stress. Strain derivatives of the Gibbs free energy G, on the other hand, 
do not yield properly defined elastic constants. 

We greatly appreciate helpful discussions with L. Stixrude. The work was supported 
by the US National Science Foundation under grants EAR-99080602 and EAR-0310139, 
and by the US Department of Energy ASCI/ASAP subcontract B341492 to Caltech 
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